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STATEMENT  PE  1'KOttUM 

The  major  ob  jec t i ves  of  this  research  are  to  develop  a general  computational 
capability  for  the  prediction  of  time-dependant  turbulent  flown  which  may 
involve  appreciable  t low  separation  and  strong  viscous- i nv i ac id  interaction 
and  to  obtain  numerical  results  that  contribute  to  the  understanding  of  the 
problem  ot  dynamic  stall.  Towards  these  objectives,  a two-equation  differential 
model  of  turbulence  has  been  studied  in  detail  and  used  in  conjunction  with 
a integro-di f terent ial  approach  for  the  numerical  solution  of  separated  flow 
problems.  A number  of  turbulent  flow  problem  have  been  solved  numerically 
and  the  results  obtained  have  been  compared  with  available  experimental  data. 

IMPORTANT  RESULTS 

The  most  significant  result  of  this  research  is  that  a versatile  and 
highly  efficient  approach  has  been  made  available  for  numerically  solving 
complex  turbulent  flow  problems.  The  cornerstone  of  this  research  project 
|l  is  the  integro-di fferent ial  approach  developed  by  the  principal  investigator 

and  his  co-workers. 

V 

The  mathematical  foundation  of  this  approach  was  established  prior  to 
the  initiation  of  this  project.  A number  of  important  attributes  of  this 
approach  was  pointed  out.  The  approach,  however,  was  only  applied  to  study 
several  low  Reynolds  number  flows  involving  simple  flow  geometries.  Under 
the  present  project,  the  approach  has  been  further  developed  and  vised  to  study 
highly  complex  separated  flow  problems  at  high  Reynolds  numbers.  The  important 
attributes  of  the  approach  has  been  conclusively  demonstrated.  User-oriented 
packages  ol  computer  codes  have  been  prepared  for -various  types  of  separated 
flows.  H 

A most  important  attribute  of  the  integro-di f terent ial  approach  is  its 
ability  to  confine  the  solution  region  to  the  vortical  region  of  the  flow. 

For  problems  where  the  tlow  Reynolds  number  is  not  small,  the  vortical  region 
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occupies  only  a very  -small  part  of  the  total  fluid  domain  and  possesses  a 
small  length  scale.  The  remainder  of  the  flow  is  irrotational  and  possesses 
a much  larger  length  scale.  A suitable  grid  system  must  have  grid  points 
closely  spaced  in  the  vortical  region  in  order  to  have  sufficient  solution 
resolution  and  accuracy.  With  conventional  numerical  methods,  the  irrotat 
ional  and  vortical  regions  are  treated  together.  It  is  difticult  to  design 
a grid  system  that  possesses  closely  spaced  grid  points  in  the  vortical  region 
and  yet  does  not  contain  an  exceedingly  large  number  of  grid  points  tv'  repre- 
sent the  entire  flow.  This  difficulty  is  one  of  the  major  motivations  for 
the  development  coordinate-transformation  and  finite-element  techniques  in 
recent  years.  With  these  techniques,  it  is  possible  to  design  an  "expanding" 
grid  system, the  spacing  between  grid  points  increasing  with  the  distance  f rom 
a solid  boundary.  In  this  manner,  the  number  of  grid  points  in  the  irrotatioual 
region,  which  is  often  far  from  the  solid  boundary,  is  minimized.  These  tech- 
niques alleviate,  but  do  not  eliminate,  the  difficulty  of  having  tv'  deal  with 
two  vastly  different  length  scales  concurrently.  In  particular,  with  increasing 
Reynolds  number,  the  length  scale  in  the  vortical  region  decreases  so  that 
this  difficulty  becomes  more  acute.  The  ability  of  the  integro-di f ferent ial 
approach  to  confine  the  solution  field  to  the  vortical  region  obviously  re- 
moves this  difficulty.  This  unique  ability  has  been  implemented  under  the 
present  project.  The  deve lopinent  of  this  approach  has  been  brought  to  a rea- 
sonable stage  of  soph i s t icat ion  and  completion  so  that  various  two-dimensional 
flow  problems  can  bo  solved  routinely  now  using  widely  accessible  computers. 
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For  example,  several  problems  involving  unsteady  laminar  flows  past  airfoils 


have  been  studied  numerically  under  separate  projects.  The  unsteady  flows  P 
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included  those  associated  with  an  airfoil  in  oscillatory  pitching  motion  and  hi 


those  associated  with  an  impulsively  started  airfoil.  Several  airfoil  config- 


urations have  been  studied,  including  many  high  angle-of-attack  cases. 
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The  details  of  these  numerical  studies  are  described  in  recent  articles 
listed  at  the  end  of  this  report  (Articles  10,  17,  19).  These  details  are 
not  repeated  here.  Representative  computing  requirements  are  stated  here, 
however,  to  indicate  the  present  state  of  development  of  the  integro-di f fer- 
ential  approach.  Typically,  for  a unsteady  laminar  flow  past  an  airfoil, 
less  than  10  minutes  of  CDC-6600  CPU  time  is  required  to  advance  the  numeri- 
cal solution  by  one  dimensionless  unit  of  physical  time,  i.e.,  the  time  inte- 
val  during  which  the  airfoil  advances  by  one  chord  length  relative  to  the 
freestream.  Based  on  the  experience  acquired  under  this  project,  it  is  con- 
cluded that  most  two-dimensional  laminar  flow  problems  now  can  be  solved  rou- 
tinely and  economically  using  the  integro-dif ferential  approach.  In  fact, 
for  these  problems,  there  is  no  critical  need  to  further  improve  the  compu- 
tational efficiency  already  achieved.  Nevertheless,  it  is  worthy  of  note 
that  several  additional  refinement  has  been  explored  for  the  purpose  of  drastic 
further  reductions  in  computing  requirements  (Article  22).  These  refinements 
are  studied  mainly  with  turbulent  flows  and  three-dimensional  flows  in  mind. 

The  most  important  refinements  studied  thus  far  are  the  separate  treatment 
of  the  boundary-layer  region  and  the  recirculating  region  (both  regions  are 
vortical)  and  the  use  of  finite  Fourier  series  rather  than  polynomials  in 
the  numerical  approximation  procedures.  It  is  expected  that  these  refine- 
ments when  fully  implemented  will  enable  the  solution  of  three-dimensional 
problems  using  modern  digital  computers  such  as  the  CDC  7600. 

In  addition  to  its  superior  solution  efficiency,  it  has  been  shown  hy 
numerical  illustrations  and  by  analysis  that  the  i.ntegro-dif ferential  approach 
offers  superior  solution  accuracy  (Articles  l,  2,  10,  15,  16,  21),  removes 
certain  major  difficulties  associated  with  the  numerical  treatment  of  boun- 
dary conditions  (Articles  8 and  ll),  and  is  useful  under  quite  general  cir- 
cumstances (Article  9).  It  is  useful  for  compressible  as  well  ns  incompres- 
sible flows,  time-dependent  as  well  as  steady  flows,  three-dimensional  as 
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well  as  two-dimensional  flows,  turbulent  flows  as  well  as  laminar  flows. 

Under  the  present  project,  considerable  progress  has  been  made  in  the  use 
of  the  integro-di f ferential  approach  for  the  solution  of  turbulent  flow  pro- 
blems . 

It  is  well  known  that  it  is  not  realistic  in  the  reasonable  future  to 
consider  the  numerical  solution  of  the  complete  time-dependent  Navier-Stokes 
equations  for  turbulent  flows,  including  all  significant-size  turbulent  ed- 
dies. Such  computations  would  require  a formidable  amount  of  computer  time 
and  data  storage.  The  development  of  a reasonably  accurate  models  of  tur- 
bulence that  simulates  the  flow  in  important  aspects  is  therefore  a more  mean- 
ingful and  rewarding  goal.  In  this  regard,  it  is  well  known  that  an  accurate 
and  general  model  of  turbulence  for  separated  flows  does  not  exist.  For  time- 
dependent  flows,  the  current  state  of  knowledge  is  especially  lacking.  Until 
recently,  the  lack  of  highly  efficient  numerical  procedures  has  severely  limit- 
ed  the  ability  of  researchers  to  test  and  calibrate  various  proposed  models 
of  turbulence.  Because  of  the  necessarily  empirical  foundation  of  turbulence 
modeling,  such  tests  and  calibrations  must  be  extensive.  The  remarkable  so- 
lution speed  and  accuracy  of  the  integro-di fferent ial  approach,  implemented 
under  this  project,  have  been  utilized  in  an  extensive  study  of  various  tur- 
bulence models  as  applied  to  separated  flow  problems. 

The  development  of  an  integro-dif ferential  formulation  for  time-dependent 
separated  flow  has  been  completed  under  this  project.  This  formulation  par- 
titions the  problem  into  its  kinetic  and  kinematic  aspects.  The  kinematic 
aspect  relates  the  mean  vorticity  distribution  at  any  instant  of  time  to  the 
mean  velocity  distribution  at  that  instant.  This  relation  is  expressed  as 
an  integral  representation  for  the  velocity  vector  and  permits  an  explicit, 
point  by  point,  computation  of  the  velocity  field.  The  kinetic  aspect  con- 
sists of  differential  transport  equations  describing  the  change  of  vorticity, 
turbulence  energy,  and  turbulence  energy  dissipation  through  diffusion,  ad- 
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vection,  production,  and  dissipation.  A turbulent  coefficient  of  diffusivity 
is  modeled  in  terms  of  the  turbulence  energy  and  the  turbulence  energy  dissi- 
pation. In  this  manner,  the  differential  two-equation  k-e  model  of  turbu- 
lence is  incorporated  into  the  integro-dif ferential  procedures  previously 
established  for  the  solution  of  laminar  flow  problems.  The  details  of  this 
new  formulation  for  turbulent  flows  are  presented  in  Article  15.  For  steady 
flows,  a standardizea  package  of  computer  code  was  prepared  (Articles  13  and 
20)  for  internal  flows  within  arbitrarily  prescribed  boundaries.  This  code 
required  only  the  input  of  grid  point  locations  and  boundary  conditions. 

No  additional  efforts  for  handling  complex  boundary  shapes  are  required. 

This  package  is  user-oriented  and  highly  efficient.  It  has  been  utilized 
in  calibrating  various  turbulence  models  in  the  steady  state  limit. 

Encouraging  results  have  been  obtained  for  a variety  of  turbulent  flow 
problems,  including  flows  associated  with  finite  plates,  cylinders,  cavities, 
and  airfoil  sections.  It  has  been  found  that  the  k-e  two-equation  differen- 
tials model  yields  numerical  results  in  good  agreement  with  available  experi- 
mental data  under  various  circumstances.  Some  of  the  results  obtained  have 
been  published  (Articles  15,  16,  and  21).  Additional  results  are  being  pre- 
pared for  publication. 

In  addition  to  the  above  summarized  efforts  in  computing  laminar  and 
turbulent  separated  flows,  a generalized  theory  for  aerodynamic  forces  and 
moments  have  been  developed  under  a separate  project  (Article  23).  This  theory 
offers  a great  deal  of  insight  into  various  aerodynamic  phenomena  and  permits 
various  mechanisms  of  production  of  lift  and  drag  forces  to  be  identified 
and  interpreted.  The  theory  has  been  utilized  in  the  study  of  unsteady  aero- 
dynamic forces  and  moments  acting  on  oscillating  airfoils. 
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